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ABSTRACT 


The classical working states of the airscrew were 
studied employing the three-dimensional flow visualization 
tunnel at the Naval Postgraduate School, Monterey, Cali- 
fornia. The various states were simulated in an attempt to 
gain insight into some of the present problems in heli- 
copter rotor design. This method was extended to include 
two particular phenomena, the critical part of the vortex 
ring state and Mangler's distribution of induced velocity 
about a rotor in the helicopter mode. Using a feasible 
distribution ascertained from the flow visualization, the 
MeOMcCe@ VelOCity Gisi«ibuution was Caliculatec for several 


points in the vortex ring state. 
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TABLE OF SYMBOLS 


Specified values are used in subsequent calculations. 


steady velocity of ascent or descent 

induced velocity at blade element dr 

radical distance of element dr from rotor axis 
rotor radius = 30 feet 

number of blades 

enora of blade (considered constant) 


geometric blade angle measured from no-lift line 
= .07 radians 


angle of attack measured from no-lift line 
moe chr onanmuy 40s 


As 
we 


“ 
Ww 


angular velocity of blades, considered constant = 
Zi radians/sec 


solidity or ratio of blade area to disk area 
BC/TR = .05 


Gemcikiy Of air 
elemental coefficient of lift 
thrust along the shaft 


elemental thrust due to B blade elements at radial 
distance r 


non-dimensional rotor radius r/R 


slope of lift curve = 3 
flevecnrough rotor; we =v ror Upwash, Vo we for 
downwash 


it 





I. INTRODUCTION 


The last decade has witnessed the emergence of the 
helicopter's increasingly important role in military and 
Civil aviation. Militarily, it has become the backbone 
of the Army's tactical doctrine. The other services have 
also found it increasingly reliable in their own particular 
operations. Civil aviation has barely scratched the surface 
of the potential applications of the helicopter. It has 
been used as a shuttle service between large urban areas 
and their major airports. Recent studies have shown that 
this service could be profitably extended to short and 


[gc Ceeretiume, in campetitimem wlth the nassender 


wEaium a 
jets. Like the proverbial tortoise and the hare, the heli- 
copter, which can transport the passenger from door-step 
to door-step, is often a faster mode of transportation than 
the jet liner. 
The demand for faster and more versatile helicopters 
has renewed interest in low speed aerodynamics after a long 
emphasis on supersonic, hypersonic, and space applications. 
The flow field of a helicopter rotor is extremely diffi- 
cult to model mathematically. The periodicity of the rotor, 
areas of reverse flow, and the non-uniformity of the inflow 
are only a few of the problems which are confronted when 


trying to calculate flow and performance parameters. 


ay 





One tool used to aid in these calculations is flow 
Visualization. From a visual study of the flow in a rotor 
system, a feasible mathematical expression for the different 
variables may be determined. In this paper flow visuali- 
zation is used to examine the working states of the rotor 
and to assist in the calculation of induced velocity in 
the vortex ring state, that is to say, recirculatory flow 


at negative speeds of advance. 








II. FLOW VISUALIZATION 


With an increasing interest in STOL/VSTOL aircraft, 
some of the unsolved low speed aerodynamic problems, which 
have fallen by the wayside in the space age, are receiving 
more attention. As discussed previously, one of the more 
useful tools in selecting a mathematical model is flow 
visualization. With this in mind a flow visualization 
facility has recently been completed at the Naval Post- 
graduate School. This facility is modeled after the one 
described in Reference 5. 


Several methods were attempted to obtain the cloud used 
as the visuelization mediwm. The m@set effective was found 
to be the vaporization of silicone-based turbine solvent. 
This fluid produces a dense white coherent cloud. Low- 
pressure air forces the fluid to drop on a heated surface 
in a locally manufactured cloud generator heated by a 
2400-watt hot plate. As the vaporized fluid collects in 
the top of the generator, it is forced by the low air 
pressure through a nozzle directed over the blade. The 
rotor is a one-bladed device driven by a variable speed 
electric motor. The speed of the blade and tunnel as well 
as control of the visualization devices is centrally located 
near the observation window. A strobe light with a variable 


flasher allowed the blade to be viewed once every revolution. 


Since the vortices are being shed at essentially the same 


frequency as the blade rotates, they also appear to remain 
stationary in the free stream. Pictures of the flow were 
taken with a Polaroid camera using ASA 10,000 film. 

Figure 2 shows a schematic of the tunnel. Figure 1 is a 
picture of the symmetrical cross-section blade mounted 
horizontally. The blade is counterbalanced to insure rotor 


vibrational stability. 





Figure 1. Rotor Used in Horizontal 
Segments of Visualization. 
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Pit. wens GUAUERT EMPIRICAL CURVE 


In 1925 C.N.H. Lock and others proposed to augment the 
"vortex theory of airscrews" with empirical data. Since 
this theory breaks down for a range of working conditions 
between zero speed of advance and some negative speed of 
advance, a two-bladed airscrew was tested at angles of 
attack below the stall. The data from this airscrew was 
used to fill the gap in the range where the theory did not 
apply. In 1926 H. Glauert extended Lock's analysis to 
cover the complete range of working conditions. The data 
which eventually evolved is commonly referred to as Gléuert's 
eel camyc. This curve is expressed in aemms af two 


non-dimensional coefficients, F and f£, defined by the 


following relationships: 


ar... ea: 

ao a Amo Oe woe (1) 
dia 2 

aa 470 es E (2) 


Me. is the speed of ascent or descent, v is the induced 
velocity at the blade element, and r is the distance of the 
element from the rotor axis. The results of these two 
investigations are given in Figure 3. 

The helicopter rotor in some instances acts as a pro- 
peller, notably in the hover and in vertical ascent or 


descent. The conventional anti-torque or tail rotor is 


10 





subjected to the same type of flow when the helicopter is 
in sideward flight. Therefore the empirical data is of 
interest in the study of helicopter aerodynamics. In the 
following pages the three working states of the rotor will 
be discussed and flow visualization pictures will be used 


towsubstantiate the discussion. 


A. THE NORMAL PROPELLER STATE 
The normal propeller state is the region on Glauert's 


empirical curve (Fig. 1) which starts at a point where 





Figure 3. Glauert's Empirical Curve. 
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1/£ = 0 and where 1/F has one of several values. According 
to Glauert's curve 1/F = 2.0; by the momentum theory 1/F = 
1.0. In practice 1/F is between these two extremes, nor- 
mally around the value of 4/3. The normal propeller state 
extends to the right as values of 1/F and 1/f increase. 

This state can be represented as shown in Figure 4, T being 
Piewenrust of the rotor and u the flow through the disk. 
iitctially, in the hover mode, is = 0 and u= v. As power is 
increased, the helicopter hacen wee a speed uae The 

flow through the rotor is now Bx + v and, from the momentum 


theory, is vy + 2v downstream. As the speed of ascent 
Votv 
Vo 


In the limiting case of u/V, = ] Sey fe o, Therefore 





increases, the value of Aad 1) since u/V, = 1+v/V.- 


f- FP=- TW — HA ana wv — W 
-— -_ _ ~~ whe el v e 


Oo 


ees 


0 ee 
Figure 4. Diagram of Normal Propeller State. 
The photographs in Figure 3 were taken with the rotor 


arranged to simulate the normal propeller state. The rotor 


axis is mounted parallel to the free stream which flows from 


r2 





left to right. The blade is rotating clockwise when viewed 
along the free stream vector. The thrust and induced 
velocity vectors are in the same direction as depicted in 
Figure 4. The tip vortices are obvious as is the vena con- 
tracta. The decrease in slipstream area downstream indicates 
an increase in velocity. In Figure 5-b the start of cyclonic 
tails which originate in the eye of the vortex are plain. 
Although not obvious in the second picture, these tails or 
trailing vortex lines extend downstream of the rotor in a 
sprial pattern, the visible series of tip vortices being 


simply the result of vortex flow around each spiral vortex 


lamer 





a. Gs 


Figure 5. Normal Propeller State. 
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1. The Hover 

As discussed previously, one of the limiting cases 
of the normal propeller state is the point where 1/f = 0 
and 1/F = 2. This is the hover mode and presents a severe 
challenge to the aerodynamicist in his attempt to predict 
rotor performance. The induced velocity is extremely 
important in this region because in a hover the induced 
power is approximately two thirds of the hover power for a 
helicopter with a conventional figure of merit. Previous 
attempts to calculate the rotor performance have not been 
very accurate because the vorticity from a blade is not 


swept away from the disk and strongly influences flow on 


the preceding blade. As the number of blades increases, 





Figure 6. The Hover. 
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this effect becomes more pronounced. Figure 6 shows the 
blade mounted horizontally, with VA = 0 to simulate the 
hover mode. In the pictures the trail of vorticity is 
quite evident and is still in close proximity to the plane 
of the disk, even though the blade is past the beam posi- 
tion. With a four-bladed rotor, the next blade would be 
strongly affected by this vorticity. The pictures also 
show the vortices extending to the periphery of the disk 


and indicates the vena contracta under the rotor. 


B. 8 2HE VORTEX RING STATE 

If the helicopter is in a hover and power is decreased, 
it begins to descend under power. This is the vortex ring 
state and is represented empirically on Glauert's curve as 
the branch extending from the hover point, which has pre- 
viously been discussed, to the so-called point of ideal 
a@eeosecrmeaon, that is, 1/f = 2.0 and 1/F = 0. Schematically 


it is represented by the following diagram: 


< 
{ 
—| 


Figure 7. Diagram of Vortex Ring State. 
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In this case the thrust vector is in the same direction 

as the relative velocity in a descent. Initially, as power 
is reduced, the induced velocity v is much larger than the 
sinking speed and u/V is negative and large. As the 
Sinking speed increases, u/V, becomes smaller and the 
induced velocity does not travel as far upstream before it 
is dissipated by the free stream. With further increase in 
Ve the tip vortex flow is driven downstream (with respect 
to a, and is then drawn into the rotor flow. This recir- 
culatory phenomenon starts at the rotor tips and moves 
inboard. At some sinking speed, ee, the induced velocity 
can not escape from the plane of the rotor disk and the 
whole disk becomes embraced in turbulent flow. 


in Figure 6 the anyie of attack 


allow the test section to simulate a rotor in the vortex 
ring state. The free stream moves from left to right and, 
this time, the thrust vector points in the same direction. 
Figure 8-a simulates the rotor at a very low sinking speed. 
The tip vortices are clearly seen. The cloudy background 
indicates that the induced velocity has been dissipated 

and the flow is downstream outboard of the rotor disk. 
Figure 8-b shows the critical part of the vortex ring state. 
The blade is enveloped in a turbulent flow, the pattern of 


which is fluctuating periodically. The phenomenon results 


in a leoccmeomusthrusct and severe Vibration of the rotor. 
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a. lane 


Figure §, Vortex Ring State. 


C. THE WINDMILL BRAKE STATE 

The last of the working states is the windmill brake 
state and occurs when the blades are extracting energy from 
the air. It is the upper branch of Glauert's curve from 
the ideal autorotation point to 1/f = 1/F = ». Now 
u= Vo - v and u/V, =o, Ae ve increases, U/ Nie iliey 


V 
fe) 


Schematically this state is shown below. 


Figure 9. Diagram of Windmill Brake State. 
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Figure 10 is an extention of the previous pair of 
photographs except that the tunnel speed is increased above 
the value which caused the blade to be in the critical part 
of the vortex ring state. Both pictures clearly show that 
tip vortices are again formed and that the area downstream 
is greater than that upstream or the converse of the stream 


tube embracing the disk in the normal propeller state. 





a. lone 


Figure 10. Windmill Brake State. 
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IV. MANGLER'S DISTRIBUTION 


Glauert's assumption that the induced velocity is 
distributed uniformly across the disk is fairly restrictive. 
In 1948 Mangler calculated an induced velocity distribution 
which was more general. He assumed that there was an infi- 
nite number of blades, allowing the rotor system to be 
replaced by a circular disk. He also assumed that the tip 
speed ratio was large. This rules out investigation of 
the region around the hover mode, which has already been 
pointed out as a critical area. Mangler selected a plaus- 


ible loading over the disk (Fig. 11) and calculated the 


Fh 


Correspenaing Giestribution o 


u. 


ke 
re) 


disk loading 
So 


non-dimensional 
a 








Oo 4 2 3 : 5 : 
on-dimensional radius 


Figure 11. Non-dimensional Disk Loading. 
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The results are shown in Reference 4 for different values 

of disk incidence. The most significant deviation from 
Glauert's original assumption is that, since the disk is 
producing littye@ere 1S circulation around the disk. Thus 
there is an upwash in the front of the disk. This phenomenon 
is shown clearly in Figure 12, there being an upwash in the 
front of the disk followed by a downwash. Notice also the 


vortices in the downwash through which the blades are 


passing. 
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Figure 12. Mangler's Distribution. 
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Vo eewecuUlinx TION OF THE LINEAR STEADY STATE VELOCITY 
IN THE VORTEX RING STATE 

With the aid of smoke flow visualization, it is now 
possible to chose a plausible mathematical model in an 
attempt to calculate the induced velocity in the vortex 
ring state. These velocities will be determined by con- 
sidering the data from Glauert's empirical curve and 
applying a method of analysis similar to that used by 
Bennett in Poneeence 1s 


Goeneacer Glauert's definition for the thrust 


Socrricvents: 
aT Z ray 
ar =" 47501 & (i) 
ar _ 2 
ae = 4mrpv. i ) 


u is the resultant flow through the elemental annulus and 
is expressed as ale in an area of upwash and v - ve i 
an area of downwash, both expressions being positive. To 
calculate the induced velocity, a new thrust coefficient 


F will be defined: 
aT = 
ater 4mrpv F (3) 


Dividing equation (1) by equation (2) yields 
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Dividing equation (1) by equation (4) gives, 


a = 2 (5) 


Z 
i: 2 
Vv 

Re-writing equation (4) and expanding u, noting that 

a _ 2 
ae VE Biv = ve) 

io. 2} 2 

ee ee 20/0 aay /V (7) 


or 


B= 1-2/£/F + £/F (8) 


Selving for 1/F 


2 
ene ae et gyall/ 5 (9) 


rae 


For the range of values for 1/f and 1/F in the vortex ring 


state, the value of 1/F can be determined. 
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Now Glauert's curve becomes three dimensional. The projec- 


tions are shown in Figure 13. The relationship between 1/F 
and 1/f is approximated by a fourth degree least squares 


polynomial, fit. 


A 3 2 
Se le + b(F) + c(F) + a(Z) +e (10) 





he 


a: XQ [ 2 
f 
Figure 13. "Three Dimensional Glauert's Curve." 
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where a= -2.96 





b = 12.96 
ec = -16.76 
ad = 9.68 
e = 12368 
2 
ab Vv 
selving for £, recaliing that == 
if 2 
V er 
O 
6 4 2 
zZ V V Vv 
i eS guinea, pt SS - 6 (11) 
3 vs y ee V EO xg 


Substituting into equation (2) and multiplying the terms in 


the brackets by a 





8 6 4 
V V V 
aE Amor c fy eho) bio Oo . avy 2 (12) 


eomevderean expression for the thrust in terms of airfoil 


and dynamic pressure considerations, 


=C BCp (wr) 7 z (13) 


g 


Consider an elemental blade section, as in Figure 14. 


Cc Aa = A(0+$) (14) 


Ns 


Vo V 
€: A(6+ pe ) (15) 





L 


Substituting equation (15) into equation (13) and equating 


the result to equation (12), recalling that 


v_ 8 V V ; 
2 _O. 
foe 


Row“rAe 





Wr 





6+ ‘3 Wee 
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Vo- V 
Figure 14. 
_ RowAe 
Let Ci = — 
oo 
C wr 6 + a 
C WRXO + CLV, - 


u°F = V oe 

Oo 
2 

Oe 
V 

Fr-2>7Ft 
Oo 

V. =vA 

Oo 





a = 
aVuv + v )F = Ve £ 


Elemental Blade 
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Section. 

Vo" cS oe Z 

V Ny. 
(18) 

6 4 

b ve (© Vo 2 

a2 4° = 2 

F V F ov 
(19) 
(20) 
(245) 
Zz) 
(23) 








where 
2FF 


a ee 
PRSLF-£F 


Le is the constant speed of ascent or descent, and v is 


the average value of induced velocity. Substituting 


equation (23) into equation (18) 


23 B23 =a 
CiF WROX+ CF vA - Cli Fv = 
v-F - a®v" - pPv7A° - cP2v7A°- aF va (24) 
Finally 
Ro 8 ae ee ” (as) 
il 2 3 
where 
Se (Fo+ 2.96A°- 12.96PA° + 18.76F2A* - 9.68F°A*) 
K,=C Flom) 
D 1 
i 2 Se ae ash: 
3 i 


Knowing the blade characteristics, the location on the 
three-dimensional thrust curve, and the non-dimensional 
distance along the blade, the induced velocity may be cal- 
culated. Some typical values are given below in tabular 


form and Guaphacally in Figure 15. 


/ £ W/E 1/F A Vaverage v2 comments 
0 ZO 2a0 0 24.5 0 hover 
5 1.0 2-591: ~42 So. 15-1 eondi tion one Fig.t5 
ays. 2.69 -47 36.23 T7e2 
nad) 250 ~26 330. LenS 
os B25 Zeo/ nos S026 ZO. 0 econdition two Figs 
2.0 0 2.0 1.0 29.7 2527) Jdeal autorotacion - 
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VI. CONCLUSIONS 


Flow visualization is an extremely useful tool which 
should be further exploited now that the facilities are 
available. More sophisticated methods of photography could 
be employed to record the vortices. Velocity measurements 
could be made in the wake with a device such as a hot-wire 
anemometer. These measurements, coupled with the photo- 
graphs, could help clarify some questions about the complex 
flow in the rotor wake, as well as other flow problems. 

Tip vortices exist in all the working states of the 
rotor just as an airplane wing has trailing vortices in all 
regimes Of riot. @al Nive trewer me@e the: vortices gan nar 
escape from the plane of the rotor before the next blade 
comes along. The flow over the blade is, therefore, 
extremely complicated and rotor performance is difficult 
EQ Calculate. 

At some point in the vortex ring state, the steady state 
velocity of descent prevents the induced vortices from 
leaving the blade. The whole rotor is enveloped in vortex 
rings which are symmetric about the axis. This condition 
is unstable and these vortices are shed in a non-uniform 
manner. 

Mangler's distribution is shown to be at least quali- 
tatively correct. There is an upwash in front of the disk 


and the blade cuts through its own vorticity closer to 
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the hub. This blade-vortex interaction could cause unde- 
sirable vibrational effects. 

The induced velocity distribution for several points in 
the vortex ring state was calculated. One of these points, 


"jdeal autorotation," is not an empirical point, as is the 
others, but a theoretical point from Reference 3. The 
values obtained appear to be optimistic, that is, half the 
blade in upwash and half in downwash. As reported in 
Reference 1 the value where Vo = v is much closer to the 
tip. This calculation takes into consideration the profile 
drag of the blade. As disk loadings become higher, the 
figure of merit approaches one, the profile drag will 


become less important and the distribution will approach 


- 2 gre) a <— ee ee | Soon a 
eae Vaart case Cepurmvea Lm Ss 
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feasible distribution ascertained from the flow visuali- 
zation, the induced velocity distribution was calculated 
for several points in the vortex ring state. 
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